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Abstract 

In the Randall- Sundrum model where the radion is stabilized by a Goldberger-Wise 
(GW) potential there is a supercooled transition from a deconfined to a confined 
phase at temperatures orders of magnitude below the typical Standard Model critical 
temperature. When the Higgs is localized at the IR brane the electroweak phase tran- 
sition is delayed and becomes a strong first-order one where the Universe expands by 
a few e-folds. This generates the possibility of having the out-of-equilibrium condition 
required by electroweak baryogenesis in the electroweak phase transition. We have 
studied numerically the region of the GW parameter space where the theory is consis- 
tent and the latter possibility is realized. We have found that in most of the parameter 
space the nucleation temperature is so low that sphalerons are totally inactive inside 
the bubbles. The condition for sphalerons to be inactive after reheating imposes an 
upper bound on the reheating temperature that is weaker for heavy Higgs bosons so 
that the out-of-equilibrium condition seems to favor heavy over light Higgses. The 
condition for sphalerons to be active outside the bubbles puts an upper bound on the 
number of e-folds at the phase transition, roughly consistent with the critical value 
required by low-scale inflation to solve the cosmological horizon problem. 



1 Introduction 

The Standard Model (SM) of electroweak interactions is being tested to high accuracy 
at present (e.g. Tevatron) and past (e.g. LEP) colliders and the physics it describes will 
be probed at the future LHC at CERN. In spite of its impressive experimental successes 
the SM has a number of theoretical and experimental drawbacks that make it difficult 
to be considered as a fundamental theory. In particular, from the experimental side 
the SM does not have a good candidate for the Dark Matter of the Universe while on 
the theoretical side it exhibits a number of problems as e.g. the unnatural fine-tuning it 
requires to accommodate the weak scale along with heavier ones, as the Planck scale. In 
order to solve some of these problems several extensions of the Standard Model have been 
proposed so far. The most popular of them, the minimal supersymmetric SM extension 
(MSSM) which solves the big hierarchy problem and provides interesting candidates to 
Dark Matter, shares some of the other SM problems. A completely different solution 
to the hierarchy problem was proposed in Ref. [1] where the Higgs field is localized in 
a four-dimensional (4D) brane inside a five-dimensional (5D) space with a warped fifth 
dimension: the warping red-shifts large (Planck-sized) scales to TeV ones, thus solving the 
big hierarchy problem. More recently, warped models with matter in the bulk have been 
widely discussed with the aim of finding a solution to the little hierarchy problem as well. 
An incomplete list includes Higgsless theories [2J, in which the electroweak symmetry is 
broken by boundary conditions and no Higgs scalar is needed, and models of Gauge-Higgs 
Unification in warped space [3]. The latter scenario is also interpreted, via AdS/CFT, as 
a calculable version of Composite-Higgs and appears particularly promising. 

An interesting feature of the SM is that it contains baryon number violation by non- 
perturbative effects (sphalerons at finite temperature) that have the capability to generate 
the baryon asymmetry of the Universe at the electroweak (EW) phase transition via the 
so-called EW baryogenesis (EWBG) mechanism [4j. Unfortunately one of the require- 
ments of EWBG, a strong first-order phase transition, was soon proved not to be fulfilled 
in the Standard Model (see e.g. [5j), thus preventing this appealing possibility. While 
the possibility of having a strong enough first-order phase transition is marginal in the 
MSSM [6] we will prove in this paper that the Randall-Sundrum (RS) theory provides in a 
natural way a supercooled first-order electroweak phase transition that can accommodate 
the mechanism of EWBG provided an extra source of CP-violation is generated. The 
main point of the mechanism is that the radion delays the electroweak phase transition 
till temperatures much lower than the electroweak one where the order parameter of the 
Higgs effective potential <j){T)/T is large. We summarize in the rest of this section the 
essential points of this mechanism. 

At finite temperature there are two stationary solutions of the 5D gravity partition 
function [7J . One is the usual RS geometry with two branes located on the UV and the IR 
points respectively. The other one corresponds to AdS-Schwarzschild (AdS-S) geometry 
where the IR brane is replaced by a black hole horizon. Using the AdS-CFT correspon- 
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dence both phases correspond in 4D respectively to confined and deconfined strongly 
interacting gauge theory. Because conformal invariance is only spontaneously broken in 
the confined phase the free thermal energy of the AdS-S phase is lower than that of the 
RS phase and thus the latter is metastable. However after explicitly breaking conformal 
invariance of the RS phase, e.g. by the introduction of a 5D field which creates a potential 
for the radion stabilizing the distance between the UV and IR branes as in the Goldberger- 
Wise (GW) model [8], the free energy of the confined and deconfined phases becomes equal 
at a given critical temperature and for lower temperatures the phase transition from the 
AdS-S to the RS phase can proceed. We have studied numerically the process of bubble 
formation (of the IR brane out of the black hole horizon) and found that Q, though the 
critical temperature lies in the sub-TeV region the actual nucleation temperature is orders 
of magnitude below, essentially due to the flatness of the radion stabilizing potential. Thus 
the radion phase transition is a supercooled one and it corresponds to a strong first-order 
phase transition. 

On the other hand the Higgs is localized on the IR brane: it can be either exactly 
localized or exponentially localized, in both cases it corresponds to a composite object 
from the holographic point of viewQ In the AdS-S phase, then, the Higgs is deconfined 
and so from an effective theory point of view it can be described as being in the symmetric 
phase = 0. In the RS phase the Higgs field appears (confines) and so the possibility of 
having <p ^ opens up. The radion supercooling prevents then the Higgs phase transition 
to proceed at typical electroweak temperatures. When the barrier between the deconfined 
and the confined phases is sufficiently low to be overcome by a quantum jump at low 
temperatures the ratio (fi(T) /T can be large and the electroweak phase transition becomes 
a strong first-order one. In the present article, for simplicity, we restrict ourselves to 
the case of an exactly localized Higgs. In the conclusion we will comment on how this 
mechanism can be applied to other scenarios like Higgsless or gauge-Higgs unification 
theories. 

The contents of this article are as follows. In section 2 we describe the two phases, RS 
and AdS-S, from the holographic point of view and provide the free energy of both phases. 
In section 3 we add the Goldberger-Wise bulk field in order to stabilize the radion potential 
and the Higgs localized in the IR brane that provides a vacuum expectation value (VEV) 
to spontaneously break the electroweak symmetry. We identify the region in the parameter 
space where there is a stable minimum fixing the radion VEV and where the back-reaction 
on the RS metric is small and the theory is consistent. In section 4 the phase transition 
is analyzed in some detail using both numerical methods and analytical approximations. 
In most of the available parameter space the Euclidean action is dominated by 0(4) 
symmetric bubbles and the nucleation temperature corresponds to a few e-folds of inflation. 
In section 5 we apply the results of the supercooled electroweak phase transition to the 

1 Some preliminary analyses of the phase transition can be found in Refs. [71 [9l 110]. 
2 The rest of the SM matter fields (fermions and gauge bosons) can be either localized in the IR brane 
or they can be propagating in the bulk of the fifth dimension. 
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required conditions for EWBG. Sphalerons inside the bubbles are inactive in all cases 
because the nucleation temperature is far lower than the SM phase transition temperature. 
Sphalerons have to be active outside the bubbles and this leads to an upper bound on the 
number of e-folds of inflation, as iV e < 26, roughly consistent with the number of e-folds 
(N e > 30) required to solve the cosmological horizon problem with TeV scale inflation [11] . 
A second condition is obtained from requiring sphalerons not to erase the generated baryon 
asymmetry after reheating, which implies a maximal reheating temperature. Counter- 
intuitively the larger the Higgs mass the weaker this constraint becomes. This latter 
condition, for a given number of e-folds of inflation, translates on upper bounds on the 
IR scale (fJ-Tev) -obtained from the minimum of the GW potential- as a function of the 
Higgs mass. This leads to values of fiTeV in the TeV range at least for a moderate number 
of e-folds of inflation. Again the larger the Higgs mass the weaker this condition. Finally 
section 6 is devoted to our conclusion and a list of open problems. 

2 Holography of the two Phases 

The conjectured AdS/CFT correspondence [12] relates extra-dimensional (stringy) models 
of gravity to AD gauge theories formulated on the boundary of the extra space. One is often 
interested in the small string length limit in which the gravity side is well described by 
effective (super-)gravity and, according to the correspondence, the AD theory has large 't 
Hooft coupling. The inverse loop expansion parameter for gravity is dual to the number TV 
of colors, so that weakly coupled extra-dimensional models are dual to large— N strongly 
coupled AD theories. In fact we will never be able to exit the strong ('t Hooft) coupling 
regime with our field-theoretical description because, unlike the more familiar QCD case, 
both the confining and deconfining phases are strongly coupled in the dual AD theory we 
are discussing. 

We will consider the Randall-Sundrum (RS) set-up [1], i.e. 5D gravity with negative 
cosmological constant and two four-dimensional boundaries which are called UV and IR 
branes. Holographically, this corresponds [13J to a AD conformal field theory (CFT) 
coupled to gravity in which the conformal symmetry is also spontaneously broken. The 
spontaneous breaking is an IR effect and its scale \i is associated to the position of the 
IR brane. The presence of the UV brane, on the contrary, explicitly breaks the conformal 
symmetry as it couples the theory to gravity. This occurs at an UV scale k = L^ 1 , 
where L is the AdS radius, which also represents an UV cut-off for the CFT. At the 
technical level, a more precise statement of the correspondence is as follows [12| 114). The 
partition function Z4 [M] of the AD gauge theory formulated on the space M with metric 
g is computed as a constrained partition function for 5D gravity 



where the subscript U G = g" means that the 5D metrics G(x,z) one integrates over are 
only those which induce a given AD metric g on the UV boundary. In Eq. (12. lj) S grav is 




(2.1) 
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the standard 5D gravity action with boundaries and will be better specified below. Let us 
now briefly remind what Eq. (|2.ip implies for the standard interpretation of the RS model. 
The l.h.s. is the partition function of the CFT calculated with a cut-off k, which contains 
not only conformal-invariant terms but also local (non-conformal) counterterms such as 
for instance an Einstein term for the AD metric g [14J. In the RS set-up the metrics 
one integrates over have no constraints at the UV, so that the RS partition function is 
obtained by integrating again the r.h.s. of Eq. (|2.ip on all the AD metrics g(x). On the 
l.h.s. this corresponds to make gravity dynamical with a AD Planck scale set by the CFT 
cut-off k. However in our context gravity is not required to be made dynamical and all 
we need is Eq. (|2.1|) in its Euclidean version. 

The thermal partition function of the AD theory is indeed obtained by considering a 
space M = S 1 x T 3 where the circle S 1 represents Euclidean time and the 3-torus T 3 
ordinary three-space. The metric is Euclidean, the length of the circle is (3 = 1/T (T is 
the temperature) and we will call V the volume of T 3 . At finite temperature, as noticed 
in [7J, the 5D path integral in Eq. (|2.ip has two stationary solutions at the classical level. 
The first one is the Euclidean version of the standard RS geometry with compactified 
time. It is 

ds 2 RS = ^ [p 2 dr 2 + V 2 / 3 d£ 2 + dz 2 ^ , (2.2) 

where the extra coordinate z is in the [zjjv = L, zm = fJ,^ 1 ] interval, r is a temporal 
variable with unit period and £ span a square 3-torus of unitary volume. The second 
solution is the so-called AdS-Schwarzschild (AdS-S) space, with metric 



dsAdS-S 



£ (ft (1 - dr 2 + V 2 '^ 2 + T -^- I ) . (2.3) 



Also in this case we chose the extra coordinate z £ [L, Zh] so that the UV brane is located 
at zjjv = L, but the end of the space now is not the IR brane, but the black hole horizon 
Zh where the space (|2.3p has, in general, a conical singularity. However, by choosing 

z h = ~Ph , (2.4) 

7T 

we can put the deficit angle to zero and, under this condition, the space is completely 
regular and it is a true solution to the Einstein equations. The time periodicity f3h in 
Eq. (|2.3p is determined by Eq. (|2.ip which states that at the UV the length of the time 
circle must be equal to j3 = T _1 . It is 



2 r 



Pl = Y [Wl + 4vr 4 (Z/r)4j . (2.5) 

The physical meaning of the two solutions is quite simple [7j. At finite temperature, 
one has two classical minima of the 5D action, and two disconnected semiclassical (loop) 
expansions could be carried out around each of them, leading of course to completely 
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different physical results. Prom the 4Z) dual point of view, the RS and AdS-S minima 
correspond to two different phases, which we identify, respectively, with the confining and 
deconfining phases of the gauge theory. It follows from Eq. (|2,ip that the free-energies 
in the two phases are proportional to the (effective) actions of the two gravity solutions. 
Depending on the temperature, and neglecting for the moment the dynamics of the phase 
transition in an expanding Universe, the physical minimum is the one with lowest 5D 
action or (equivalently) 4D free-energy. 

At the classical level one computes the free-energies by plugging the solutions (12. 2|) 
and (|2.3[) into the action. The latter reads 



= -4M i 

^ grav 



(R + I2k 2 ) + / 2K UV 
M JdMuv 



+ f (2K IR + 6k) 

JdM IR 



(2.6) 



where R is the 5D curvature scalar, M ~ Mp the 5D Planck mass and -fQyy(/m the 
extrinsic curvatures of the two boundaries dMuy^jny Following the conventions of |15j . 
in Gaussian normal coordinates 

K = 1 ah K ah = X - 1 ab d r)lah , (2.7) 

where r\ is the coordinate normal to the boundary and 7 the induced metric. The equations 
of motion associated to the action (|2.6p are the 5D Einstein equations in the bulk plus 
the Israel junction condition 

K^ = -k^, (2.8) 

at the IR brane only. One has no junction condition at the UV because the UV metric is 
fixed in (|2.ip and should not be varied when working out the equations of motions from 
the action. 

While both spaces (|2.2|) and (|2.3p solve the bulk Einstein equations only the RS ge- 
ometry solves the Israel junction condition (|2,8p . thanks to the fine-tuned choice we did 
in Eq. (|2.6p for the IR brane tension. On the contrary the AdS-S metric (|2.3p does not 
solve Eq. (|2.8p . meaning that an IR brane cannot be put in that geometry. Thus strictly 
speaking both spaces (|2.2p and (|2,3p are not solutions of the same gravity action. In fact 
for the AdS-S to be a solution the IR boundary terms (last terms of Eq. (|2.6p ) should be 
dropped. One might think however that the IR brane is hidden by the black hole hori- 
zon in the AdS-S space so that the IR terms of the action do not effectively contribute. 
More precisely one should think of the IR brane in RS as a remnant of a more compli- 
cated solitonic-like geometry which solves the equations of motions of a more fundamental 
underlying theory of gravity. The same theory should also allow for a second solution, 
without the brane, which would correspond to the AdS-S space. 
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The free-energy of the two phases is readily computed. By plugging Eq. (|2.2p into 
Eq. (|2.6p we get for the confining (RS) phase 



F RS 



PV 



S 



grav 



[RS] 



3 lA 



-2A(ML) s k 



(2.9) 



which is just the usual RS bulk cosmological constant which we could have canceled by 
adding a suitable UV brane tension term to Eq. (|2.6|) . For the deconfined phase we must 
remove the IR terms from Eq. (|2.6|) and plug Eq. (|2.3|) . One gets 



Fa<is-S 



1 

pv 



~~mr Sgrav [AdS S] 



3, A 



-24A; 4 - 4vr 4 T 4 ( 1 + O 



8(MLyk 
7r 4 r 4 



A: 4 



Ph P 
(ML) 3 , 



(2.10) 



where an expansion for irT <C k has been performed. 

As a general rule of AdS/CFT, the loop expansion on the gravity side corresponds to 
a large number of colors (large-AQ expansion in the gauge theory. The squared coupling 
constant for 5D gravity is 1/(ML) 3 and naive dimensional analysis suggests that the 
expansion parameter is given by the squared coupling divided by a loop factor. The 
precise AdS/CFT relation 1 



1 

iV2 



(ML\ 
16vr 2 



(2.11) 



perfectly matches this interpretation. Let us now compare the above equation with the 
free-energies (|2.9|) and (|2.1U|) . The RS phase is confining and indeed -Frs only contains 
a vacuum energy term to which all the N 2 (— 1) confined gluons can contribute. The 
vacuum energy scale is fixed by the UV cut-off k of the AD theory. The AdS-S phase, 
on the contrary, is deconfining and the N 2 (— 1) massless gluons are present in the plasma 
everyone giving its standard ~ T 4 contribution. The vacuum energy, which is the same 
in both phases, can be always shifted to zero by adding an UV brane tension term to the 
action (12.61). 



3 Adding the Golberger-Wise and Higgs fields 

Equations (|2,9p and (|2.10p immediately uncover a problem: at any non-zero temperature 
FAdS-S < Fjis so that the RS space is metastable and no phase transition can ever 
occur from the AdS-S to the RS phase. This is because the conformal symmetry is only 
spontaneously broken in the RS phase and it is completely unbroken in the AdS-S one. 
Being spontaneous, the conformal breaking indeed leads to a massless Goldstone boson, the 
radion, which arises in 5D as a modulus of the RS geometry corresponding to the distance 

3 The relation (|2.1ip can be changed by numerical factors depending on the precise theory. We will use 
Eq. (|2.11[) as a definition. 
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among the two branes [13]. Having fixed the location of the UV brane, the radion field 
controls the position fi^ 1 of the IR brane. Given that the radion is a Goldstone, it cannot 
have a potential and this is why it does not contribute to the vacuum energy and then to 
the free-energy (|2.9p at leading order in the large-iV expansion. Since the radion potential 
is flat and fi undetermined, the temperature is the only dimensionful quantity and there is 
no other scale the temperature can be compared with (the cut-off k is not relevant here). 
Stated in another way there is no scale which could fix the temperature of transition T c . In 
fact starting in the AdS-S phase and lowering the temperature the system would remain in 
that phase forever. The same occurs for the first holographic deconfining phase transition 
discussed by Witten [16] in the case of unbroken conformal symmetry. 



3.1 The Golberger Wise field 

Independently of the above considerations an explicit breaking of conformal symmetry 
is required anyhow to avoid the presence of an exactly massless scalar. The so-called 
Goldberger-Wise (GW) mechanism [17] permits to stabilize the distance among the branes 
and give a mass (and a potential) to the radion by simply adding a 5D scalar field The 
Euclidean action is 

Sow = [ [G MN d M *d N $ + m 2 $ 2 ] + / (£/*(*) + A: 4 <m) , (3.1) 
JM JdMiR 

where ST\ amounts to a correction to the IR brane tension in Eq. (|2.6[) . We do not need 
to specify the IR Lagrangian Cm since it only plays the role (as in [1?]) of enforcing the 
IR boundary condition 

${x,z IR ) = k^VL (3.2) 

Holographically the GW field corresponds [13] to an explicit UV breaking of the con- 
formal symmetry due to an operator O of conformal dimension d = A + e where 

2 r2 



m 2 L 



e = y/A + m 2 L 2 - 2 ~ — ^— , (3.3) 

and we have assumed m 2 <C k 2 to perform the expansion. In AdS space the squared mass 
of a scalar can be positive or negative and then e can have both signs. To have a solution 
of the hierarchy problem |e| must be small (< 1/10) and we are going to assume this in 
the following. Adding the GW term (|3.ip to the action modifies the correspondence of 
Eq. (|2.ip in a very simple way. On the l.h.s. one has the partition function of a different 
AD theory which gets modified by the introduction of the term 

A $(x) O , (3.4) 



£e+3/2 



to the AD Lagrangian. In the above equation A is a dimensionless coupling, k is the UV 
cut-off and <3?(x) a non-dynamical AD source field. On the r.h.s. of Eq. (|2.ip we now have 
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the new 5D field to integrate over. As for the metric one has to consider a constrained 
path integral. The value of $ on the UV brane is 



$(x,zuv) = $ = k 3/2 



VQ 



(3.5) 



where we fixed the source $ to a constant value and then converted the operator (|3.4p 
into a deformation of the CFT. One could also have made dynamical, and added a 
Lagrangian for $ at the UV as in [T7j. The role of those terms, as for Cir, would simply 
have been to enforce the boundary condition (|3.5p and hence, in practice, to convert $(x) 
back into a non-dynamical field. The two approaches are then completely equivalent. 

To find the background configuration in the presence of the GW scalar one has to 
solve the coupled Einstein and Klein-Gordon equations which come from the total action 
Sgrav + ^gw anc ^ impose the boundary conditions (|3.2|) and (|3.5|) . This is in general 
a non-trivial task and exact analytic solutions can only be obtained in some particular 
cases [IS]. The problem greatly simplifies if we can neglect the backreaction of the GW 
scalar on the 5D metric. In this case one simply has to solve the Klein-Gordon equation 
on the corresponding unperturbed gravitational background, provided by Eqs. (|2.2p and 
(|2.3p . and this provides an approximate solution of the original coupled equations. In the 
RS case the actual value of /i (i.e. the position /i" 1 of the IR brane) is the one which 
minimizes the value of the action on the corresponding solution. By plugging the solution 
into the action one gets, in addition to the vacuum energy term of Eq. (|2.9|) . a potential 
Vgw(p) f° r the radion whose minimum fixes the value of the (stabilized) modulus. The 
result is [IT] 

V GW (^) = ev 2 k 4 + /j, 4 [(4 + 2e)(«i - v (fi/k) e ) 2 - evj + 5T X ] , (3.6) 

in the very good approximation of neglecting (/x/£;) 4 <C 1. 

For the potential (|3.6p to have a (non-trivial) global minimum one must require 



6T ± 

tor e > , — 2~ < 6 i 



ST, 



for e < 0, -(4 + e) < —± < e 



(3.7) 



If the previous conditions are satisfied the global minimum of Vgw(/*) is located at 



1A 



UTeV = k [ — | , where X 

>0, 



4 + e + sign(e) /e(4 + e) - 4^ 



iSTi 



4 + 2e 

It is also useful to have a formula for the value of the potential at its minimum. It is 



(3i 



V{HTev) = tvlk A + fJ-TeV v l 



2 + e 



— - sign(e) Je(4 + e) - 4— 



(3. 
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Equation (|3.8[) shows how the GW mechanism provides a solution to the Hierarchy 
Problem. Since k is associated to the Planck scale and UTeV to the Electroweak (TeV) 
one, an enormous fine-tuning is expected "a priori" in order to get [iTeV out of microscopic 
parameters of order k. However by making a sensible choice of |e| and vi/vq not too far 
from 1 Eq. (|3,8p can easily give rise to the required 16 orders of magnitude among fiTeV 
and k. Solving the Hierarchy Problem is the greatest success of this class of models, and 
considering too high values of |e| would spoil this success. If we at most allow for a two 
orders of magnitude hierarchy among v\ and vq we get the upper bound: 

kl < § - (3.10) 
For |e| = 1/20, which is the reference value we will mostly use in the following, vi/vq ~ 

10 si S n ( € ) 4 / 5 . 

For the above manipulations to make sense the GW field must provide a small back- 
reaction on the metric. To ensure this condition, the contribution of the GW field to the 
energy-momentum tensor must be compared with the one coming from the bulk cosmo- 
logical constant. This leads to the constraints 

tt 2 v 2 3 1 //A 2c 

L < / 



iV 2 |e(4 + e)| + e 2 X 2 \k 

TT 2 vf 3 

~W < | e (4 + e )|+ e (4 + e)-4^' (3,U) 

where the |e(4 + e)| term in the denominators comes from the fact that, depending on 
the sign of e, the stronger constraint arises from the T^ u or T zz components of the 5D 
energy-momentum tensor Tmn- The small back-reaction condition also implies another 
constraint, which was considered in [17] and ignored in [9]. It comes from comparing the 
RS brane tension in Eq. ()2.6f) with the 5T\ term of Eq. (13. lj) . Remembering that the latter 
has been treated as a correction, we must impose 

|<5Ti| < 2A(MLf = -AriV 2 . (3.12) 

Depending on the sign of e significantly different constraints on v\/N arise. When 
e is positive (^/k) 2e is small and the first bound in Eq. (13. lip is the relevant one. The 
two bounds could only be comparable for very small values of e. In the following we will 
however mainly be interested in negative e. In this case the first bound in Eq. (|3.1ip is 
weak and the second bound becomes the relevant one. As for the latter v \ drops from the 
second line of Eq. (|3.1ip one gets, noticing that 5T\ < 

IflTil < 12(ML) 3 = -\n 2 , (3.13) 

which is a factor of 2 stronger than Eq. (|3.12p . Considering Eq. (13. 7p . Eq. (I3.13P can be 
converted into a bound on v i as 

7T 2 V 2 3 



< — . (3.14) 



N 2 " 4|e 
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The above condition is always stronger than the one coming from the first line of Eq. (|3. 1 1 [) . 
Indeed, (n/k) 2e > 1 and it is possible to show, using Eq. (|3.7|) . that X is smaller than one. 
We can then conclude that for e < Eq. f)3. 13j) is the only small-backreaction condition, 
although one should however account for the consistency conditions in Eq. ([37 
Focusing on the negative e case it is convenient to define the variables 



4tt 2 |(5Ti| Air 2 \e\v\ 
3 N 2 ' U ~3 N 2 



(3.15) 



The perturbativity constraint (|3.13j) and the consistency conditions (|3.7p can respectively 
be expressed as 

6><i, v < e < 4 7j e ^ - ( 3 - 16 ) 

|e| 

As far as e is small, the last inequality is substantially irrelevant and all our constraints 
can be expressed in terms of two parameters only. 

In order to discuss how the GW field affects the deconfined (AdS-S) phase, following 
[7], we have to construct an "off-shell" version of the AdS-S geometry (j2.3|) by defining a 
"Hawking temperature" = l/(jrzh) and relaxing the regularity condition (|2,4p . A 4D 
field Th(x) is associated to the parameter T/j. Very much like the radion field /jl(x) controls 
the position n~ l of the IR brane in the RS phase, Th(x) controls the position Zh of the 
AdS-S black hole horizon. Differently from the radion, however, T/, has a potential even 
in the absence of the GW field. Having detuned the condition (|2.4p . indeed, the space 
develops a conical singularity which, once regularized with a suitable spherical cap, gives 
rise to a potential for with a minimum at T h = T for T <C k). Taking also 

into account the GW contribution, the potential [i.e. the deconfined (AdS-S) phase 
free-energy] is 

F d (T h ,T) =E + ^fTt - nWTTl + ^fuX 2 {^-^ , (3.17) 
where we subtracted the vacuum energy term 

#0 = -VGW^Tev) , (3.18) 

in order to cancel the cosmological constant term [Eq. (|3.6p at the minimum] of the confined 
(RS) phase. At small |e| the above formulas reduce to 



F d (T h ,T)~E + tt 2 N 2 

with 



\(i + v -)t£-ttI 



(3.19) 



3N 



2 



E o = J^2 V(9v-" 2 )\e\l4ev ■ (3- 2 °) 
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Notice that the above expansion holds for \e\ 1 ^ 2 -C 1 and fails for 0>i/. In this approxi- 
mation Fd has a minimum at 

T mm = T_ /g 21) 

h 1 + I//2' V ; 

and the free-energy at the minimum is 

^ 2 iV 2 T 4 

Fr»(r) = Bo -— IrTw . (3.22) 

The introduction of the Th(x) field associated to the horizon position also permits 
to develop an intuitive understanding of how tunneling among the two vacua occurs [7]: 
starting in the AdS-S phase at = T™ m the horizon starts to move away from the UV 
brane since reaches 0. At this point the horizon has disappeared and the metric ()2.3p 
reduces to the one of the infinite AdS§ space, which coincides with the RS space (|2.2[) for 
fi = 0. The IR brane now appears from infinity and /i grows from zero to its equilibrium 
value UTeV- The radion, in this schematization, is the only field of the RS phase which 
undergoes a variation during the tunneling process and it is only possible to justify this 
approximation (at a not too high temperature compared with UTev) if the RS KK gravitons 
are significantly heavier than the radion. In this case it is energetically expensive for the 
KK fields to move and the transition can be studied in an effective theory for the radion 
alone. With the aim of checking that indeed m r < tukk ~ ^^TeV i n & U our allowed 
parameter space, and for completeness, we report here the radion Lagrangian [19|, [T3l [8] 

3N 2 1 3iV 2 1 

Cr = 2^2 (a ^ )2 " Vgw{ ^ = 2^2 [ ( ^ )2 " + " 'J ' (3 ' 23) 



where 



m 2 = 2(2 + e)v (AX + AeX - 2eX 2 - AX 2 ) n TeV 

~ %y/-e{0-v)v& eV . (3.24) 

Level curves of m r /rriKK in the region of parameters allowed by Eq. A3. 16H are shown 
in Fig. [H and we see that the radion is always sufficiently light for our approach to the 
transition to be meaningful. 

Introducing the SM fields, as we will do in the following, does not dramatically change 
the picture of the transition, we will however have to include a second light scalar, the Higgs 
field, on the RS side. The transition becomes, in principle, a problem of tunneling in a two- 
dimensional field space, which is numerically much more complicated than the standard 
one-dimensional one. We will argue, however, the Higgs contribution to the bounce to 
be negligible and studying the one-dimensional radion tunneling will be sufficient for our 
purposes. 

We conclude this section by observing that small-backreaction is not at all necessary 
for the GW mechanism to work, as shown in [18J. As long as |e| is sufficiently small 
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0.2 0.4 0.6 0.8 1 

6 



Figure 1: Level curves of m r /rriKK equals to 0.25 (red short-dashed line), 0.19 (purple dashed), 
0.13 (black dot-dashed) and 0.06 ( green solid) in the plane (9, is). We have fixed e = —1/20. 

the GW mechanism will still provide a solution to the Hierarchy Problem even in the 
presence of a significant backreaction. From this point of view the small backreaction 
conditions we impose only have a technical motivation, as they ensure the consistency of 
our manipulations. All calculations done in this section could in principle be repeated 
taking backreaction into account. One should consider the coupled Einstein and Klein- 
Gordon equations, find the two solutions and substitute them back into the action. This 
would lead to exact expressions for the free energies of the two phases at leading order in N. 
For what concerns the dynamics of the phase transition, on the contrary small backreaction 
is more than a technical assumption. The point is that, as shown in Fig. [I] this condition 
ensures the radion to be light. If this were not the previously discussed, the 

schematization of the transition by the radion tunneling would be meaningless and one 
would need to go to the true extra-dimensional gravitational instanton which connects 
the two vacua. This is most likely beyond the reach of our field-theoretical description of 
the two phases, as one should face the problem of describing in some way the process of 
the IR brane hitting the horizon and disappearing into the black hole. 

3.2 The Higgs field 

Let us now introduce the Higgs field in our theory. As in the original RS scenario we take 
the Higgs doublet h to be localized at the IR brane. The (Euclidean) action is then 



4 We are assuming here for simplicity that there is no Higgs-radion mixing from a term like 



SM 




(3.25) 




)h ] h, i.e. that £ ~ 0. 
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where the Higgs doublet h has a tree-level potential 

Vh(M) = A (^ _ ^ 2 ) • ( 3 - 26 ) 

For the rest of SM fields we will not consider a particular scenario although we will assume 
that gauge bosons and (at least part of) the SM fermions[§ are propagating in the bulk. 
We will also assume that for low enough energies (where all Kaluza-Klein excitations are 
decoupled) the zero modes are effectively described by the SM. 
Using now the RS metric (|2.2p one can write the action as 



/<M(; 



S#m= I d A x{ (JL_Y\ Dh f + f_E_Y Fo(/i) „) + ...) (3.27) 

\/J>TeV J \HTeVj 



where the Higgs field is redefined as 

i, ^-h. r.. _/'^- 



, v ~ 246 GeV . (3.28) 
k 

Notice that Higgs fields and mass parameters (or vacuum expectation values), are red- 
shifted with a power of the factor fiTev/k according to their dimension. 

The IR brane terms only contribute to the RS phase and hence, in the present model, 
the Higgs field is not present at all in AdS-S. The holographic reason for this is that the 
Higgs particles are entirely composite objects which cannot be present in the deconfined 
high-temperature plasma. 

Light SM fields, on the contrary, are most likely fundamental particles. Including their 
contribution the AdS-S free energy becomes 



TT 2 N 2 T 4 



7T 



2 



F d = E - — : ; XQ - — g*, T 4 , (3.29) 
4 (l + iz/2) 3 90 yd ' v ; 



while in the RS phase we have 



7T 2 ^4 



F c (f,, cf> c , T) = VrsQi, <Pc, T)+E Q - —g* c T , (3.30) 

where g* ~ 100 {g* d ) is the effective number of light SM degrees of freedom in the RS 
(AdS-S) phase, 4> c = y2h° and Vrs{^j 4>c T) indicates the radion-Higgs potential. At the 
tree-level it reads 

\ 4 \ 



Vrs(h, <t>c) = Vawfji) + I - J — -(<f>i-v 2 y. (3.31) 

\HTeV J 4 

As we can see from Eq. (]3.31|) the typical size of the SM potential is Vsm ~ v 4 while that 
of the GW potential is much larger, i.e. Vow ~ Mtw Therefore the total potential 

Vrs = V GW [1 + 0{v A /v A Tev)] (3-32) 



5 Some fermions could also be localized in the IR brane as the Higgs fields. 
6 As well as other particles possibly localized on the IR brane. 
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can be approximated, as far as the phase transition from AdS-S to RS is concerned, by 
the GW radion potential. In fact at the typical radion transition temperatures the SM 
potential is a little perturbation of the radion potential and should not alter the bounce 
solution and the corresponding Euclidean action that governs the phase transition. The 
radion potential then produces a supercooling of the system and the electroweak phase 
transition takes place at much lower temperatures than the typical EW ones and makes it 
to be much stronger. We postpone the numerical analysis of this behaviour to section 5. 

In the approximation of Eq. (13.311) the free energy is minimized for [i = fiTeV and 
d) r = v , and at the minimum 



TT- 



—g* T 
90 yc 



(3.33) 



Although irrelevant for what concerns the radion phase transition, one-loop corrections to 
the Higgs potential are important at high enough temperatures to determine the VEV of 
the Higgs field. For [i = fiTeV the one-loop Higgs potential is just the standard one since 
we have assumed a SM-like low energy dynamics, and is given in the following. 

In the MS-scheme the one-loop Coleman- Weinberg potential depends on the renormal- 
ization scale Q through the ratio m 2 /Q 2 while the thermal corrections similarly depend 
on the temperature through the ratio m 2 /T 2 . We can then write the SM potential at 
/i = UTeV as 



VSM 



64tt 2 



n-i[mi{ 

--W,Z,h, x ,t 



Q 2 



Ci 



+ 



jt4 
2^ 



y~] niJ B [m 2 /T 2 
-W,Z,h,X 



+ n t J F [m 2 t /T 2 



(3.34) 



where the thermal fermionic (bosonic) function Jpm) are given by 

1 ± e ~^ x2+y2 



Jf(b)(v 



dx x log 



(3.35) 



the constants C% derive from the MS renormalization scheme we have used, and are the 
degrees of freedom of the % field, that is 



Cw = Cz 



6 



Ch — 



a 



nw = 6, nz = 3, n h = 1, n x = 3, n t 



2 

-12 . 



(3.36) 



The only fermion that contributes sizeably to the effective potential for <j) c ^ is the 
top quark t since the other fermions are very weakly coupled to the Higgs doublet. For 
practical purposes we fix Q = mt{v) which is nearby the Higgs VEV. As for the thermal 
corrections it could have been possible to use a better approximation including plasma 
effects [20], however we think that for the purposes of the present paper the one-loop 
approximation is good enough. 
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4 The Phase Transition 



At high temperatures (early times) a deconfined (strongly-coupled) plasma fills the Uni- 
verse, which is described by the Friedmann-Robertson- Walker metric ds 2 = dt 2 — a 2 (t)dx 2 
and expands according to the Friedmann equation 

a „ 8tt 

a= H =\l3M 2Pd > (41) 

where pd is the energy density of the plasma and can be easily extracted from the free- 
energy in Eq. (|3.29p . As the Universe expands, it cools down and a confining phase 
transition can occur when the free-energy of the deconfined phase equals the confined 
one. This happens at a critical temperature T c which is obtained by equating Eqs. (13.290 
and (13331 : 

T c 4 = ^{6v- V 2 )\e 

where Ag* = g* — g d and g d is the number of light degrees of freedom in the AdS-S phase 
that equals, in the model we are considering, to the Higgsless SM degrees of freedom. At 
small enough iV and large enough u, can become negative. This is due to the fact that 
the T term in Eq. (|3.33j) becomes larger than the one in Eq. (|3.29j) and in this case, given 
that Eq is positive, the deconfined phase is never stable, the Universe stays forever in the 
confined one and no phase transition occurs. However for N > 3, and using Ag* = 4 in 
our model, this situation is avoided provided that the constraints (|3.16p are fulfilled 0. In 
most of the allowed parameter space the critical temperature is a bit smaller than fiTeV'- 

4.1 Completion of the transition 

Formally, the phase transition begins as soon as the temperature drops below T c . It will 
then be convenient, to fix ideas, to choose the origin of time as T(0) = T c and rescale 
the space coordinates so that a(0) = 1. There is a simple one-to-one relation among 
temperature and time, T(t) = T c /a(t) and, using Eq. (|4.ip 

f = —TH(T) = -2LyJ^ Pi (T). (4.3) 

For T <T C there is a non-zero probability for bubbles of the confined phase to form, the 
rate of bubble nucleation for unit physical volume can be expressed as 

A(T) = A(T) e" 5(T) , (4.4) 

7 If some of the fermions, like for instance the top quark, are localized to the IR, Ag* would be larger. 
However, all our results are substantially insensitive to its precise value. 



1 2Ag* 
(l + v/2) 3 ~ 45iV 2 



IJ'TeV 
7T 



(4.2) 
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where the exponential suppression is due to semiclassical tunneling and the bounce action 
S(T) will be computed in the following. The prefactor A can be estimated on dimensional 
grounds, A ~ T a T^ with a + b = 4. In particular for very low temperatures, as those 
that will be found in this work, a = 0, b = 4. The nucleation rate is dominated by the 
exponential so that this rough estimate of the prefactor will be sufficient for our purposes. 

To have an idea of how fast the transition proceeds we follow [21] and write down the 
probability that one point of space remains in the old (deconfined) phase at the time t. 
Assuming the bubble expands at the speed of light 



p = exp 



exp 



47T 

~3~ 



dtia 3 (ti)A(T(ti)) 



o 



47T 

T 



T J l 




(4.5) 



where we used Eq. ()4.3p to change variable in the time integrals. The nucleation rate, 
due to the exponential factor in Eq. (|4.4p . undergoes exponential variation in the course 
of time, i.e. temperature, so that there could be regions of temperature in which the dT\ 
integral in Eq. (|4.5p only receives negligible contributions. Using dimensional analysis (we 
are considering now T ~ T c ~ TeV) to estimate p<i and the dT^ integral in Eq. (|4.4p (which 
we could as well compute using Eq. (|3.29p without changing the result which follows) we 
find that the region of temperatures for which 

4 

« 1 for T ~ 1 TeV 



S(T) > B = 4 log 



M P 



PTeV 



148 



A(T) (M P 



rp4 



(4.6) 



cannot contribute appreciably to the integral if T is not orders of magnitude smaller than 
pTeV ~ TeV. 



We will soon see that, for T ~ T c , the condition (14. 6p is always fulfilled in the case 
we are considering, so that the phase transition does not effectively start at T = T c . As 
discussed above, it is very unlikely to find a point in the new phase, the Universe is filled by 
a supercooled deconfined plasma and keeps expanding with the energy density associated 
to Eq. (|3.29p . It is worth remarking that, as the temperature drops below 



V(^-^ 2 )|e|(l + v/2f 



3 ffJ>TeV 



TT 



(4.7) 



the cosmological constant starts to dominate over radiation and an inflationary epoch 
begins. The total number of e-folds of inflation is ~ log Ti/T n where T n , to be defined 
below, is the temperature at which the transition effectively occurs. In the inflationary 
epoch we can use p<i = Eq in Eq. (|4.5p and write down more explicitly the probability of 
transition as 

4tt f Ti dTt A(Ti) ( M P 1 



p(T) = exp 



El 



T 



(4- 
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where we neglected the [Tj,T c ] part of the integral. Once again we can check that tem- 
peratures for which the condition (|4.6p holds do not contribute to the integral. In our 
model, the bounce action S(T) monotonically decreases with temperature so that if (a) 
S(T — > 0) > B the action never crosses the critical value while if (b) S(T — » 0) < B we 
can define the nucleation temperature T n as 

S(T n ) = B . (4.9) 

The case (a) corresponds to the "slow" phase transition considered in [22] and, despite of 
the fact that the probability (|4.8p will eventually go to zero, bubbles will not percolate, the 
phase transition will never end and we will have an empty and cold Universe. Avoiding 
(a) will give us a strong constraint on the allowed parameter space of the model and, in 
particular, on the maximum value of the number of colors N. 

Let us now focus on case (b). For T\ = the nucleation rate is always too small for 
the transition to start but for T\ = T n the argument of the integral in Eq. (|4.8j) becomes 
of order one. It is important to remark that the bounce action has a power-like behaviour 
with temperature so that the rate grows exponentially. When the temperature crosses 
T n the argument of the integral suddenly passes from being much smaller to much bigger 
than one. The integral suddenly becomes large, the probability goes to zero and the 
phase transition ends. During the time the phase transition lasts, therefore, the Universe 
undergoes a negligible cooling, i.e. a negligible expansion. On the scale of the evolution of 
the Universe the phase transition is instantaneous: bubbles expand, collide and percolate 
as if the Universe were static. 

We can illustrate this behaviour with a simple toy model. Since the integral (14. 8|) 
is dominated by temperatures very close to T n one can assume that the action in (14. 4|) 
behaves linearly with T, i.e. S{T/T n ) = a + (3T/T n where a + P ~ B is the condition (j4.9j) 
and the slope j3 has been computed numerically and yields values (3 ~ 50 — 100 (depending 
on the values of the GW parameters). We can then write for the probability (j4.8f) the 
expression 



p(T) = exp 



4tt r T * <n\ eBe _ {a+pTi/Tn) / _ t_ 

3 Jt T\ \ T\ 



(4.10) 



The behaviour of the probability p(T) is shown in the left panel of Fig. [2j We see that, in 
this particular case the phase transition happens between T = 0.88 T n and T = 0.84 T n , a 
negligible cooling that corresponds to the little expansion of 0.05 e-folds during the phase 
transition. The model considered in Ref. [22] corresponds to the case where (3 = and 
Eq. P~l"0j) reads 



p(T) = exp 



■f 60 ( Ne ~ y + 3e ' Ne " r~ 2Ne + r~ 3Ne 



(4.11) 



where eo = ex.p(B — a) and we can see that a large number of e-folds N e (T) = log(Tj/T) is 
related to the value of eo by N e ~ 3/4ireo which would imply eo -C 1. On the other hand 



18 







0.82 0.84 0.86 0.88 0.9 . 92 

T/T n 



-4 



-3 -2 -1 

logxoCzyro 



Figure 2: Left panel: Plot of the probability (|4. 10[) as a function of T/T n for T n = KT 3 !^ and 
a slope (3 — 100. Right panel: Comparison of the previous plot [red (solid) line] as a function 
of log 10 T/Ti with the probability (|4.1ip for eo = 1 [blue (dot-dashed) line] and eo = .01 [green 
(dashed) line]. 

for eo ~ 1 we obtain p < 1 after an 0(1) number of e-folds of inflation. This does not 
happen in our case since nucleation is triggered by the temperature T n where condition 
(|4.9p holds. The difference between both situations is illustrated in the right panel of 
Fig. ® 

An analogous situation was discussed in Ref. [23] and it was shown that two conditions 
must be satisfied by 7 = X/H 4 to complete the phase transition in an (inflationary) first- 
order phase transition, i) 7 > 9/4tt is required for percolation; ii) 7 ~ for almost all 
the inflationary period, otherwise bubbles are formed, grow during inflation and distort 
the CMB radiation. We have checked that in our model both conditions are satisfied at 
least for moderate number of e-folds of inflation, i.e. for N e < 6. Certainly condition 
i) is satisfied since 7 crosses from values larger to values smaller than one around the 
nucleation temperature. Condition ii) is also satisfied since H is a constant during the 
inflationary period while T varies exponentially with the temperature so that 7 varies 
very quickly near the nucleation temperature. Therefore big bubbles do not have the 
chance to develop. We did not make the numerical analysis for the very small values of 
the nucleation temperature corresponding to the number of e-folds of inflation consistent 
with the required amount of cosmological inflation (iV e ~ 30) because inflation was not 
the main issue in this paper. We will come back to this issue in a different work. 

Summarizing, we have seen that, provided Eq. (|4,9p is verified at some non-zero tem- 
perature, the phase transition completes, the new phase plasma thermalizes and the ex- 
pansion continues with the free-energy given by Eq. (|3.30p . During the transition, of 
course, the energy (or what is the same, given that the Universe is effectively static, the 
energy density) must be conserved so that the Universe will end up in the confined phase 
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with a reheating temperature Tr given by 

^9* C T R = E . (4.12) 

To obtain the above equation we equated P c (Tr) obtained from Eq. (|3.30p with Pd{T n ) 
and neglected the nucleation temperature T n <C Tj. 

4.2 Bubble nucleation: approximate analytical results 

We will now describe in detail the tunneling process from the deconfined (AdS-S) to 
the confined (RS) phase as described by the respective free-energies given by Eqs. (|3,17p 
and (|3.30p . At the critical temperature T c the two minima are degenerate and thus the 
Euclidean action S(T C ) — > oo. Below T c the first bubbles that can be formed are those 
with 0(3) symmetry, thin walls and very large radii. In this case the probability of 
bubble formation is governed by S = S^{T)/T where S3 is the Euclidean action with 0(3) 
symmetry. In the thin wall approximation the action can be given an analytic expression 
as S3 = IQtt Sf / 3(AV) 2 where Si is the surface tension that can be evaluated in the limit 
T — > T c and AV is the depth of the true vacuum relative to the false one |241I25|. Using the 
analytic approximation we will prove now that the radion phase transition will not allow 
the formation of thin wall bubbles. For that we will momentarily forget the Higgs field 
since its presence will not alter the following conclusions. A straightforward application 
of the thin wall formula for the deconfinement / confinement radion phase transition leads 
to 



S 3 (T) 64tt n 2 TeV {3N 2 \ 3/2 p TeV /T 



3 j2E~ Q \2tt 2 J (1-T 4 /T C 4 ) 2 



(4.13) 



where Eq is the vacuum energy defined in (|3.20p . Eq. (|4.13p shows that S%{T)/T goes to 
infinity in both limits T — > and T — > T c . The function S^{T)/T has in fact a minimum 
(that corresponds to the maximal probability of nucleation) at T* = T c j y/3. Using now 
the relation (|4.2h between pxeV an d T c one can write the final expression for the action 
(l4X3j) at T = T* as 



S 3 (T*) 50N 2 



T* ~ [(8 - z>|e|] 3 / 8 



1 2Ag* 



(1 + I//2) 3 45iV 2 



1/4 50N2 



where the last inequality follows from the fact that the term raised to the power 1/4 is 
an order one number and from the perturbativity bounds 9 < 1, v < 1. For instance for 
|e| ~ 1/10, S3/T does not satisfy condition (|4.9|) for N > 2 while iV > 3 is needed for a 
reliable 1/N (semiclassical) expansion of our theory. We can observe that if the values of 
the parameters (6, u) are small, as it corresponds to the perturbativity region, the values 
of the Euclidean action S3/T can be huge. 

Of course thin wall bubbles with symmetry 0(4) have Euclidean actions larger than 
those just described with 0(3) symmetry. On the other hand since there is supercooling 
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we will found that in most of the available parameter space the least action corresponds 
to 0(4) symmetric solutions. 

For 0(4) symmetric bubbles we have worked out the thick wall approximation that 
we will illustrate now and where we have also included the Higgs potential. In order to 
do that we will work out the thick wall approximation for a set of field with (in general) 
non-canonical kinetic terms as follows. Let us consider a theory with a number of scalar 
fields = (0i, (f>2, ■ ■ ■ ) with a potential V(4>) and general non-canonical kinetic functions 
/i(0). The euclidean action for 0(4) configurations can be written as 

S 4 = 2tt 2 dpp 3 ^^(0)1^ + F(0)j , (4.15) 

where we are using the notation 0' = dep/dp. The potential V has a false vacuum (say at 
= 0) where the energy is normalized to be zero V(4>) = and a true minimum at = 0o- 
The bounce solution for 0(4) symmetric bubbles is obtained by solving the equations of 
motion (EOM) 

fi<t>" + ~fi4>'i + ^jZrr^'i = IT + oxf"^' 2 no summation on i , (4.16) 
p J d<pj d(f>i 2 d(f>j J 

with the boundary conditions 4>'(p = 0) = and 0(p — > oo) = <p. 

We can now consider a bubble of true vacuum and radius R inside an exterior of false 
vacuum. Since the potential outside the bubble is zero and the 0-profile constant, only 
the region inside the bubble will contribute to the action (|4.15p . If (f>* = 0(0) is the value 
of at which the bounce starts^, 5(f) = 0* — is the total variation of inside the bubble 
wall 51? We can then approximate the action (|4.15p by 

S 4 ^ vr 2 ^ £ fi 2 SR + ^VR* , (4.17) 

where V < is a suitable average of the potential inside the bubble. We will take 
V = V((f)*) and will use the notation V(0*) = V* and so on. Notice that in this way we 
will underestimate the bounce action. 

If the bubble is thick we take 5R = R and determine the radius R c of the bubble by 
extremizing Eq. (|4.17p with respect to R. This gives for the critical radius 

V d)* 2 f* 

Rl = , (4-18) 

and for the critical action 

8 Notice that (f>* does not coincide with (f>o. 

9 The bubble wall 5R is defined as the region p where <f> varies. 
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We will apply the previous expressions to the case of the system of two fields, radion 
and Higgs fields, (//, <p c ) with non-canonical kinetic functions 



3N 2 1 



kin 



2vr 2 2 



^TeV 



2(d u nf + -^\d v hf , 



and potential 



where 



E 



f^TeV ' 



V GW {H) = (4 + 2e) 



(4.20) 



(4.21) 



+ <5Ti - ev 2 



(4.22) 



The value of the critical action is then given by 



S A 



7T 

T 



37V 2 
2vr 2 



+ 



i* 2 



-V* 



V* 



(4.23) 



which depends on the values of fi* and </>* that in turn do not necessarily coincide with 
their vacuum values fJ-TeV and v respectively. In fact the precise values of n* = fj,(0) and 
§* c = C (O) require of a numerical calculation of the bounce solution. Here we just want 
to give a plausibility argument for the phase transition to happen. 

The phase transition should happen when S4 ~ B as we have seen in (14. 9p . The actual 
value of S4 in Eq. (|4.23p depends, for fixed values of the GW potential parameters, on 
the fields fj,* and </>*. We can minimize the action S4 with respect to those fields in order 
to obtain the most favorable configuration for the bounce solution. Minimization with 
respect to fi* yields the value 



k ( — 

vo 



(4.24) 



for which the thick-wall action is found to be 



S 4 



IT' 

T 



37V 2 
2^ 2 



L*2 



+ 



^TeV 



evf 



V* 



(4.25) 



The action (j4.25H has a minimum at 0* = for light Higgs masses. However the point 
(fj,*, 0) is a saddle point of the potential and the corresponding configuration is classically 
stable. This situation is avoided for values of 0* 7^ 0. Since the function £4(0*) is 
monotonically increasing between (f>* = and 4>* = v a conservative assumption would be 



10 For the purpose of this section it is enough to just consider the tree level SM potential with radiatively 
corrected parameters. 
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to fix S^v) < B and therefore this condition will be satisfied by all configurations with 
< (f>* < v. The Si(v) < B condition can be expressed as 

3JV 2 / 2» 2 v 2 \ 2 3N 2 ( 0.4^ 2 



where the second term inside the squared is a very tiny correction reflecting the fact that 
the SM potential is a small perturbation of the GW potential. 

4.3 Bubble nucleation: numerical results 

Generally speaking, the approximation methods used in the previous section are barely 
reliable. The main limitation is that they do not provide a reliable determination of the 
nucleation temperature as it is not possible to know, a priori, which shape (thick or thin) 
the bounce solution has at a given temperature. A numerical calculation of the bounce is 
therefore needed, and in the case of a single scalar field the solution can be easily found 
iteratively by the technique of over-shoots and under-shoots 0. Clearly this can be done 
at zero and finite temperature. 

In our case we have three fields (Th, fj, and h) involved in the bounce but, as we have 
discussed before, the Higgs does contribute to the bubble action with a tiny amount. More- 
over the two remaining fields follow a single path in the two-dimensional field configuration 
space [7j and then we can consider the transition for a single field defined as 

f - ( NJ3/2/it) T h for $ < 

* ={ ( K r- \ } (4-27) 
1 ( N^j2/n\ fj, for $ > , 

which has a canonical kinetic term and the potential 1*^1 

V($, T) = F d U^/2j?,/N t) 0(-$) + F c U^lJilN $, 0, t) 8($) , (4.28) 



where is the Heaviside step function. 

Figure [3] provides an example of the numerical results we find for the bounce action 
[S = S4 or S = S3/T] with 0(3) and 0(4) symmetry in two typical points of our parameter 
space of N = 3. The comparison with analytical (thick and thin) formulas shows up a 
considerable discrepancy. The nucleation rate, at a given temperature, is dominated by the 
bounce of minimal action and, depending on the point in the parameter space, bounces 
with either symmetries can be dominant. In all cases, however, the nucleation temperature 
is quite small. For the cases with a moderate number of e-folds (N e < 6) that we have 



11 We have explicitly checked our bounce program for the potentials of Ref. [25] 
12 We are supposing for Th the kinetic term ^ri(d^Th) 2 similar to the radion one. 
This potential is a simplification of that considered in the numerical analysis. Since for T / the 

potential is discontinuous at $ = 0, it has been regularized and checked that the result is not sensitive to 

the regularization. 



23 



0.001 0.01 0.1 0.001 0.01 0.1 

T 7 UTeV T/ UTeV 

Figure 3: Solid [dashed] lines is the plot of the 0(4) [0(3)] Euclidean action as a function of Tf/iTeV 
for N = 3, e = -1/20 and, for the left [right] panel, (0, v) = (0.35, 0.1) [(0, v) = (0.6, 0.25)]. On the 
left [right] the transition occurs via 0(4) [0(3)] bubbles. The analytical thin wall 0(3) (|4.14[) is 
plotted in double-dotted-dashed and the 0(4) analytical asymptotic value (14.251) in double-dashed- 
dotted. 





Figure 4: Plot of the allowed regions in the plane (0, v) for N = 3 (upper-left panel), N = 4 
(upper-right panel), N — 5 (lower-left panel) and N — 6 (lower-right panel), with e = —1/20 and 
vq as a function of fiTeV '■ The consistency bound < v is marked by a double-dotted-dashed line. 
The region where the nucleation is allowed is on the right of the dotted line, obtained by numerical 
0(4) bubbles at T n = 0. The same upper bound, but calculated using analytical approximations, is 
the double-dashed-dotted line. The points where the nucleation occurs at some fixed temperature 
are indicated by colored solid lines. Finally, several reheating temperature curves are drawn in 
dashed. 
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analyzed numerically we have checked that the slope of the Euclidean action S(T/T n ) is 
(3 ~ 50 — 100 as it was anticipated in section 4.1. 

Most of our numerical results are summarized in Fig.d]which shows, for e = —1/20 and 
different values of N, the triangle in the 6-v plane allowed by consistency (|3,16p and, inside 
that, the analytic bound for the nucleation to happen as obtained by Eq. (|4.26p . Contour 
plots are shown for fixed nucleation temperature T n = 5-10 p^TeVi T n = 5 • 10 6 HTeV 
and T n = 5 • 10~ 2 //TeV> corresponding to an inflation of N e ~ 6, 4, 2, respectively. Regions 
of higher nucleation temperature, i.e. smaller inflation, are to the right of such level 
curves. As we can see in the allowed region a supercooling is necessary before permitting 
nucleation, and the larger N the lower T n as predicted by the analytical approximations. 
In the allowed area, 0(4) bubbles dominate over 0(3), except in the case N = 3 when 
T n > 10 -2 UTeV- At T = the GW potential, for e < 0, has a local minimum at \i = and 
zero-temperature nucleation occurs as a tunneling from this secondary minimum to the 
global one. By studying this zero-temperature tunnelinj^l we have obtained numerically 
the region in which nucleation is allowed (dotted line in the Fig. U]). We see that the 
analytical formula provides a reasonable approximation. Also in Fig. 2] level contours for 
fixed values of the reheating temperature in units of fiTeV (dashed lines) are shown. 

Until now we have considered 0(4) solutions also at finite temperature, but strictly 
speaking they have sense only if the diameter of the bubble 2R is smaller than T" 1 , the size 
of the compactified dimension at the nucleation temperature. Under that circumstance 
0(4) is a good approximate symmetry [25]. To show that explicitly we plot in Fig. [5] the 
nucleation temperature T n (without considering 0(3) solutions) and the corresponding 
Tn m = l/(2R)fiTeV as a function of 6 along the bisectrix [u = 9/2 — 1/20] of the allowed 

10° 

io J 
Iff 2 

> 
- 

H 

io J 
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'"0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 

e 

Figure 5: Lower curve: Plot of T n //iTeV as a function of 9 along the bisectrix [v = 6/2 — 1/20] 
of the allowed region in the upper- left plot of Fig. [4] (N — 3). Upper curve: Similar plot for 
Turn = l/(2R)[iTeV where R is the radius of the 0(4) bubble in the thick wall approximation. 

14 See |26| for an early study of this transition. 
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region of Fig. 1 for iV = 3B The restriction of this analysis to the bisectrix in N = 3 is 
justified by the fact that along this path the maximal nucleation temperatures are reached 
and therefore it should be the worse region for the 0(4) approximate symmetry. From the 
plot we can conclude aposteriori that 0(4) actions considered in Fig. [J] make sense. 

5 Conditions for electroweak baryogenesis 

As we have seen the electroweak phase transition in the RS model, triggered by the radion 
phase transition, is a strong (supercooled) first-order one. This is unlike the pure Standard 
Model where the electroweak phase transition is very weakly first order at any perturbative 
level 0. Since the weakness of the electroweak phase transition in the Standard Model 
was one of the main obstacles to the mechanism of electroweak baryogenesis we can see 
that the situation is completely different in the Randall-Sundrum setup. 

The first condition to be imposed for a successful EWBG is the non-erasure of the 
generated baryon asymmetry at the nucleation temperature inside the bubbles. As we can 
easily check this does not really impose any additional condition on the model parameters. 
In fact because the phase transition is strongly first order, and in particular (0 c (T n )} /T n 3> 
1, sphalerons are totally inactive inside the bubbles and the baryon asymmetry generated 
by some extra source of CP- violation (on top of the standard CKM phase) will not be 
erased in the broken phase. 

The second condition is that sphalerons have to be active outside the bubbles and in 
order to check this we need an estimate of the sphalerons rate in the symmetric (decon- 
fined) phase. The SM gauge bosons and the fermions arise from bulk fields in our scenario, 
so that they correspond holographically to fundamental AD fields which are mixed with 
the strongly coupled CFT. In the deconfined phase, then, we just have the SM fields 
(except for the Higgs) mixed with the new sector. Neglecting the mixing, which is ex- 
pected to give a small New Physics correction, the sphaleron transitions are just like in the 
Higgsless SM. In the symmetric phase (see e.g. [27]) the presence of the Higgs is irrelevant 
for the calculation of the sphaleron rate and we can then simply use the SM results. On 
dimensional grounds [28j , the sphaleron rate can be estimated as r sp /j ~ kq^T 4 where the 
coefficient k has been evaluated numerically in [27] with the result 0.1 < k < 1 and aw is 
the weak gauge coupling constant. The rate of baryon number non-conserving processes 
V B (T) is related to T sph by [29] 

Vb(T) = f N f^ > 

where Nf = 3 is the number of generations. We have to compare now the rate (|5.ip with 
15 R is calculated using the 0(4) thick wall approximation. 

16 In fact non-perturbative analyses have shown that it is so weak that it disappears (it becomes a 
continuous crossover) when non-perturbative effects are taken into account. 
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the expansion rate of the Universe given by 



X = 




(5.2) 



This gives a lower bound on the nucleation temperature provided by 



■n 



> 3 x HT 12 iV [(0i/- u 2 )\e\] 



i/i 



(5.3) 



UTeV 



which corresponds to an upper bound on the number of e-folds of inflation as N e < 26. The 
nucleation temperature widely satisfies the bound (|5.3p in the region we have analyzed. 

As we have seen the electroweak phase transition in our model is accompanied by a 
period of exponential expansion corresponding to a (few) number of e-folds of inflation. 
After the phase transition the Universe is reheated to a given temperature Tr. If the 
temperature Tr is higher than the temperature at which (p c (T)/T = 1 sphalerons will be 
activated inside the bubbles and the generated baryon asymmetry will be erased. So a 
third condition is that the reheating temperature be lower than the temperature at which 
4> C (T)/T = 1 for a fixed value of the Higgs mass. 

In order to understand the temperature at which <j> c (T)/T = 1 for the potential (|3.34p . 
we plot the value of the ratio (4> C (T))/T versus the temperature T (Fig. [6]). Here the 
three lines are calculated for a Higgs of 115 GeV (solid line), 165 GeV (dashed) and 225 
GeV (short dashed). In correspondence to these masses the critical temperature T^ M 
has been found, resulting respectively 146, 195 and 250 GeV, at which (4> C (T C )} /T c results 
approximately 0.2, 0.1 and 0.060 A remarkable feature of the SM potential is that, inside 
a temperature interval of 1 GeV, the origin passes from being a minimum to a maximum 
and during this change {4> C (T)) runs incredibly fast. Clearly this behaviour of (0 C (T)) 
remains also at a bit lower temperatures and, in fact, as we can better see in the right 
panel of Fig. [6l at temperatures of T ~ 135, 165 and 190 GeV, respectively, we address 



We can now compare the maximal reheating temperature as provided by Fig. [6] with 
the level contours of fixed reheating temperature provided in the various panels of Fig. 01 
For a fixed value of the Higgs mass we will first read from the right panel of Fig. [6] the 
maximal reheating temperature T^ ax . Second in the corresponding panel of Fig. H]we will 
look for the contour line corresponding to the same value of the reheating temperature 
for a fixed value of fJ-TeV- Then the available region in the (#, u) plane will be that to the 
left and below it. For instance if we consider mu = 140 GeV we can see from Fig. [6] that 
T rnax ^ 15Q Ge y Then if we fix ^ TeV = 1 Te y and for the A = 3 case we see in Fig. 1 

that there is a contour line corresponding to Tr = 0.15 fiTeV = 150 GeV and we can easily 

17 Here we define the critical temperature T c s M as the temperature where the non-trivial minimum of 
Vsm and the origin are degenerate. 

18 Of course for such small values of <f> c our approximation fails as a consequence of the IR problem of 
thermal field theories. 



(MT))/T ~ 1. 
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Figure 6: Left panel: Plot of (<j> c (T))/T as a function of the temperature in GeV for a Higgs mass 
of 115 GeV (solid line), 165 GeV (dashed) and 225 GeV (short dashed). Right panel: Plot of the 
temperature at which the SM potential is minimized at (cj> c (T)} /T — 1 as a function of the Higgs 
mass in GeV. 



identify the available region. Now if we increase the Higgs mass we will obtain from Fig. 
a larger value of T^ ax and therefore going back to Fig. H]the region below and on the left 
of the corresponding level curve will be bigger than the previous one. For instance if we 
consider now the Higgs mass ran = 180 GeV we can easily see that T^ ax ~ 170 GeV and 
the corresponding available region will be that below and to the left of the contour line 
Tr = 0.17 UTeV, and this region is bigger than the previous one. So for a fixed value of 
UTeV the heavier the Higgs boson the bigger the available region in the parameter space. 

Another point of view of the same condition is if we instead fix the minimal nucleation 
temperature, i.e. the maximum number of permitted inflation, in which case the maximal 
reheating temperature T^ ax given by Fig. [6] translates into an upper bound on the value 
of UTeV as can be seen in Fig. [7] where the case of T n ~ 500 MeV, equivalent to a number 
of e-folds -/V e ~ 6, has been considered. 

We can conclude that with an additional source of CP-violation the Standard Model 
with a Higgs boson localized on a brane of a warped space can generate electroweak 
baryogenesis. A condition for this to happen, as we have seen, is that the universe spans 

19 This extra source requires by itself extension of the Standard Model. Some mechanisms have been 
proposed, e.g. that in Ref. [30 , where the strength of CP violation induced by the CKM phase varies 
with the temperature. From the low energy point of view another interesting possibility would be the 
appearance of new CP-violating operators. For instance, as discussed in Ref. [31| . dimension-six operators 
~ 2 i — \H\ 2 FF might generate the observed baryon-to-entropy ratio. Finally an obvious possibility 
is to enlarge the SM particle content in such a way as to have new sources of CP-violation. An example 
is provided by the MSSM where one can have an independent source of CP- violation e.g. in the phase 
of the Higgsino mass in which case charginos (and neutralinos) can generate the baryon asymmetry. In 
our Randall-Sundrum scenario it would be enough to enlarge the SM with the corresponding fermions 
(charginos and neutralinos) as in Split Supersymmetry. However analyzing those possible sources of CP- 
violation is out of the scope of the present paper and moreover it would deserve an independent (dedicated) 
investigation. 
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Figure 7: Upper bound on \iTeV for different values of TV as a function of the Higgs mass from the 
condition that {(j) c {T R )) /T R > 1. 

less than 26 e-folds in the supercooling previous to the fast phase transition. This number 
is close to the number of e-folds required in low-scale inflation to solve the cosmological 
horizon problem [11] 

N e > 30 + log . (5.4) 

UTeV 

In our case condition (|5.4p translates into N e > 27.7 which is barely consistent with 
condition (|5.3p . We can even stress that in some regions of the available parameter space 
it might be possible to span the number of e-folds required by Eq. (|5.4p although we have 
not made a dedicated numerical analysis in this paper. 

6 Conclusion and Outlook 

In this paper we have studied the supercooled electroweak phase transition triggered 
by the radion in the Randall-Sundrum model at finite temperature in the presence of 
a Goldberger-Wise potential stabilizing the radion field and deduced the conditions for 
electroweak baryogenesis coming from the sphaleron interaction rates. We have concen- 
trated in the region of the space of parameters where the theory remains perturbative 
and the Randall-Sundrum metric is not perturbed by the Goldberger-Wise field. We have 
estimated from analytical approximations that the contribution of the Higgs field to the 
Euclidean action is negligible (at least for small Higgs masses) as compared to that of the 
radion and have worked in the approximation of neglecting the former. Our results can 
be summarized as follows: 

• In the region of the GW parameters consistent with perturbativity there is a wide 
sub-region where the phase transition is completed provided that the parameter N 
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is not too large. In fact we have found that for N > 6 it is very difficult the phase 
transition to be in agreement with perturbativity constraints. 

• In the available region there is always a certain amount of inflation which can go 
from a few e- folds to larger values. We have explicitly performed our numerical 
analysis for a moderate maximum value of e-folds (N e < 6) although in some tuned 
regions of the parameter space we might have the number of e-folds required to solve 
the cosmological horizon problem in low-scale inflationary models (N e ~ 30). 

• In parallel with the previous issue we obtain values of the nucleation temperatures 
much below typical electroweak temperatures o This makes sphalerons to decou- 
ple inside the bubbles, as required by the condition of not erasing any previously 
generated baryon asymmetry, independently of the detailed value of the nucleation 
temperature. 

• In the analyzed region bubbles percolate and there is no dangerous production of 
large bubbles that could distort CMB radiation. 

• After the phase transition the system reheats to temperatures Tn/^iTeV ~ 0.1 — 0.2. 
In the presence of supercooling the reheating temperature is essentially insensitive 
to the precise value of the nucleation temperature. 

• At the reheating temperature the condition for the sphaleron to remain decoupled 
translates into a lower bound on the Higgs mass. Therefore large values of the Higgs 
mass are favored by electroweak baryogenesis unlike in most models. 

• Sphalerons outside the bubbles should be active in order to generate baryon asym- 
metry in the bubbles wall. The condition for sphalerons to be in thermal equilibrium 
outside the bubbles translates into an upper bound on the number of e-folds of in- 
flation as A e > 26. 

• In the numerical analysis we have considered the case where only the Higgs field is 
localized on the IR brane. We have checked that similar results hold for the cases 
where also (part of) the SM fermions are equally localized which would basically 
translate into a modification of the effective number of light degrees of freedom in 
the AdS-S phase, as it can be seen from Eq. (|3.29|) . 

There is a number of open questions and further studies that are worth mentioning 
here: 

• Our numerical analysis has been based on neglecting the Higgs effective potential as 
compared to the GW radion potential. We believe that this approximation is fully 
justified for small values of the Higgs mass. However as electroweak baryogenesis in 

20 As much as we can compare our nucleation temperatures are lower than those obtained in Ref. [9]. 
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these models seems to favor large over small Higgs masses it should be worthwhile 
to evaluate the effect of the Higgs potential on the Euclidean action for such heavy 
Higgses. 

• For heavy Higgses the compatibility of the radion and Higgs interactions with elec- 
troweak precision tests is an issue and should be considered in detail [32J. 

• As the CKM phase in the Standard Model is known [33] to be insufficient for a suc- 
cessful electroweak baryogenesis it should be essential to find out some extra source 
of CP violation in the considered Randall-Sundrum model or extensions thereof. 

• As we already mentioned we analyzed numerically cases where the number of e-folds 
of inflation is moderate, N e < 6. However in some regions of our space of parameters 
we might produce a much stronger inflation and in particular the number of e-folds 
required to solve the cosmological horizon problem in low-scale inflation. It should be 
interesting to explore in detail such regions since there the radion might be identified 
with the inflaton field. 

• We have considered a simplified model where the Higgs is fully localized in the IR 
brane. In other models we could have similar effects. For instance in gauge-Higgs 
unification models the Higgs field, which arises from the extra-dimensional compo- 
nent of a gauge field, is exponentially localized at the IR brane and the situation 
seems quite similar to the one we have studied. In Higgsless models the gauge bo- 
son masses are provided by the IR boundary conditions and the only sensible way 
in which the electroweak symmetry could be restored at high temperatures is in 
the deconfined phase we have described, in which the IR brane disappears. In all 
of these models the EW phase transition should be similar to that studied in the 
present paper although dedicated studies should be worthwhile. 
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